Finite-size effects in tunneling between parallel quantum wires 
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We present theoretical calculations and experimental measurements which reveal finite-size ef- 
fects in the tunneling between two parallel quantum wires, fabricated at the cleaved edge of a 
GaAs/AlGaAs bilayer heterostructure. Observed oscillations in the differential conductance, as a 
function of bias voltage and applied magnetic field, provide direct information on the shape of the 
confining potential. Superimposed modulations indicate the existence of two distinct excitation 
velocities, as expected from spin-charge separation. 
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One-dimensional (ID) electronic systems are a very 
fertile ground for studying the physics of interacting 
many-body systems. In one dimension, the elementary 
excitations are collective spin and charge modes, the 
spectrum of which is strongly influenced by the Coulomb 
interaction Q|. An electron entering such a system must, 
therefore, decompose into the corresponding eigenmodes, 
resulting in a suppression of the tunneling density of 
states. This suppression was detected in a variety of ex- 
periments, such as tunneling from metal contacts into 
carbon nanotubes Q and resonant tunneling in one di- 
mension H . A unique feature of interacting electrons in 
one dimension, described by Luttinger-liquid (LL) theory 
0, is the decoupling of the spin and charge degrees of 
freedom, each of which propagates with a different veloc- 
ity determined by the Coulomb interaction. To date, di- 
rect experimental verification of this phenomenon is lack- 
ing. Moreover, issues pertaining to the decoherence and 
relaxation of the elementary excitations of the LL remain 
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FIG. 1: Schematics of the circuit (A) and the model (B). A 
wire of length L runs parallel to a semi-infinite wire. Bound- 
aries of the wires are formed by potentials Ui(x) and U2(x). 
Also drawn in (B) is a one-electron wave function ip{x) in 
the upper- wire confinement Ui(x). The energy and momen- 
tum of the tunneling electrons are governed by voltage V and 
magnetic field B. 



a challenge both theoretically and experimentally. 

Here we report a detailed experimental investigation 
and theoretical explanation of a set of interference pat- 
terns in the nonlinear tunneling conductance between two 
parallel wires that were first reported in Ref . mL A sketch 
of the tunneling geometry in shown in Fig7Ty(A). The 
interference appears because the tunneling process is co- 
herent to a very high degree and is due to the finite length 
of the tunnel junction. A wealth of information can be 
extracted from the interference: The pattern itself en- 
codes microscopic details of the potentials in the wires, 
while the structure of its envelope reflects the presence 
of two distinct excitation velocities per electron mode in 
the data, as expected from spin-charge separation. The 
decay of interference may also yield information on de- 
coherence processes of the elementary excitations in ID 
systems. 

Fabrication of ID quantum wires of exceptional qual- 
ity has been achieved by cleaved edge overgrowth in 
GaAs/AlGaAs heterostructures ||, and much progress 
has been made in experimental investigations of their 
transport properties j|, |(| . Recently, a method was devel- 
oped to measure tunneling between two parallel wires of 
this type j|] : The wires are formed at the cleaved edge of 
a wafer containing two parallel quantum wells, only one of 
which is occupied by a two-dimensional (2D) electron gas 
(2DEG) (cf. Ref. f§ for details on the sample). A volt- 
age bias V between the wires forces electrons to tunnel 
through a narrow AlGaAs barrier separating them [see 
Fig. |l|( A)] . Measurements of the differential conductance 
G at 0.25 K are made with standard lock-in techniques, 
as a function of V and B, a magnetic field applied per- 
pendicular to the plane containing the wires. This allows 
to determine the complete dispersion relations of the ele- 
mentary excitations in the quantum wires [Q . Prominent 
features of the measured G (V, B) can be understood in 
terms of a model which considers two infinite parallel 
wires, and accounts for the electron-electron interactions 
by means of LL theory fl |I| fjj. On the other hand, an 
observed oscillation pattern results from the finite length 
of the upper wire (UW). We show that key features of the 
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oscillations can be understood assuming that the ends of 
the UW are defined by a soft confining potential, rather 
than assuming sharp, square- well confinement. Most of 
our discussion will employ a model of noninteracting elec- 
trons, which explains the most prominent features of the 
interference patterns. Interaction effects will be discussed 
at the end. Fig. |l|(B) schematically shows the potentials 
Ui(x) and U^ipo) felt by electrons in the wires. Electrons 
in the UW are confined to a region of finite length by top 
gates at both ends of the junction. Electrons in the lower 
wire (LW) are reflected at the left end, but can pass un- 
der the right-hand gate, rendering it semi-infinite. The 
effective tunneling region is determined by the length of 
the UW, which is approximately the region \x\ < L/2 in 
Fig. |l|(B). The magnetic field gives a momentum boost 
fiQ = eBd to electrons tunneling from the UW to the 
LW, — e being the electron charge and d-the distance be- 
tween the centers of the wires; V > favors tunneling of 
electrons from the LW to the UW. 

As observed in Ref . , the differential conductance has 
strong peaks along a set of curves in the V-B plane, where 
tunneling can occur with conservation of energy and mo- 
mentum between a Fermi point in one wire and an elec- 
tron mode in the other wire. Here we focus on the region 
near the intersection of the dispersion peaks in G (V, B) 
at zero magnetic field Q and voltage, V — (Ep 2 —Ep 1 ) / e, 
necessary to overcome the difference in Fermi energies 
Ep 1 and Ep 2 of the UW and LW, respectively. (In the 
experiments, several electron modes may be occupied in 
the wires. Here we consider only the mode in each wire 
with the largest Fermi momentum along the wire.) The 
most spectacular manifestation of the breaking of trans- 
lational invariance is the appearance of a regular pattern 
of oscillations away from the dispersion curves. Fig. ^ 
shows typical examples of the patterns measured at low 
magnetic field. The lines that correspond to the disper- 
sions appear as pronounced peaks that extend diagonally 
across the figure. Additionally, we observe numerous sec- 
ondary peaks running parallel to the dispersions. These 
side lobes are asymmetric: They always appear to the 
right of the principal dispersion peaks. The result is 
a checkerboard of oscillations in region I, stripes in re- 
gion II, and no regular pattern in region III (see Fig. ||). 
When the lithographic length L is increased from 2 /jm 
[cf. Fig. |(B)] to 6 urn [cf. Fig. |(C)], the frequency in 
V and B increases by a factor of w 3. The period is ap- 
proximately related to L, d, and the Fermi velocity v-p by 
|AU|L/«f = \AB\Ld — 2-Kh/e. Upon close examination 
of Figs. |(B) and|(C), one can discern a modulation of 
the interference that causes a series of faint streaks par- 
allel to the B-axis, where the oscillations are suppressed. 
The strength of this modulation is shown in Fig. |^(A). As 
explained below, this is a moire pattern created by two 
superimposed sets of interference, each resulting from a 
distinct velocity that is present in the data. 

We base our theoretical analysis on a phenomenologi- 



cal tunneling Hamiltonian 

H tU n = Tj2 [ dx*l 1 (x)V s2 (x)e- iQx +H.c, (1) 

where ^ S i(x) is the spin-s electron field operator for the 
ith wire (i = 1 denotes the UW and i = 2 the LW) . Since 
the Zeeman energy in GaAs is small, we ignore the spin 
degrees of freedom, and characterize the electrons in the 
UW by a discrete set of energy eigenstates ip n (x). The 
eigenstates in the lower wire form a continuum, which we 
write as (fk 2 (x) , indexed by wave vector k 2 . Treating tun- 
neling to the lowest nonvanishing order in perturbation, 
we find for the current 

/cxsgn(U)|T| 2 ^|A/(7i,Q,U)| 2 , (2) 

m 

where n = n-p + sgn(U)m, with rip being the index of 
the state ip n just below Ep 1 ; J2 m ^ s a sum over integers 
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FIG. 2: Oscillations of G {V, B) at low field from 2 /im (B) 
and 6 /im (C) junctions. (A smoothed background has been 
subtracted to emphasize the oscillations.) The arrows in (C) 
mark the minima of the power P of the oscillations in mag- 
netic field as a function of the voltage V . (A) shows P for 
(B) (solid line) and (C) (dashed), for which the abscissa is 
2.3 x V. The brightest lines in (B) and (C), corresponding 
to tunneling between the lowest modes, break the V-B plain 
into regions I, II, and III. Additional bright lines in II arise 
from other ID channels in the wires and are disregarded in 
our theoretical analysis. 
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m with [sgn(U) - l]/2 < m < e\V\L/nhv F ; and M is 
the tunneling matrix element between state tp n and state 
(fk 2 , the energy of which is lower by eV. Specifically, M 
is given by 

M = J dxr n (x)e- lQx Vk2 (x). (3) 

If e\V\ is not too large, we can linearize about the Fermi 
wave vectors k Fi , and fc 2 is then given by [k^ — k F2 )v F2 = 
sgn(V)v Fl nnr/L — eV/%. 

As a starting point, we consider infinite square- well 
confinement in the region |a;| < L/2 of the UW, so that 
i/j n oc sinn7r(cc/L— 1/2). We also assume that the poten- 
tial in the LW is infinite for x < xq, so that the states in 
the LW have the form ifk 2 oc sinfetx — xq). In line with 
the experiments, we assume that k Fl 2 differ slightly by 
A&f = k Fl — k-p 2 , but neglect the difference in the uf's of 
the two wires, which is only a few percent Q . As a result, 
\M\ 2 is independent of the index n and is given by \M\ 2 = 
|M(+)| 2 + |M(-)| 2 , where M(±) =sin[K( ± )L/2]/[ K ( ± )L/2] 
and k^' = Ak F + eV j %v F ± Q . (In the limit L — > oo, the 
two terms |AfW| 2 become <5-f unctions.) Similar results 
were obtained by Boese et al. || who considered tunnel- 
ing between two infinite noninteracting wires through a 
window of finite length L j|] . 

Differentiating Eq. (g) to obtain G (V, B), we find that 
whenever the applied voltage matches a discrete energy 
level of the UW with the Fermi level of the LW, the sum 
in Eq. (Q) exhibits a step, yielding a series of (^-function 
peaks in the conductance. These peaks are not seen in the 
experiments, possibly because of smearing due to finite 
temperature and/or 1D-2D scattering in the UW. We 
may therefore disregard the discreteness of the sum and 
write I oc V \M(Q, V)\ 2 , which gives a pronounced oscil- 
latory contribution cx Vd\M \ 2 /dV to G. (When L — > oo, 
this oscillation pattern disappears and the conductance 
becomes G oc \M (Q, V)\ 2 resulting in bright peaks along 
the dispersion curves.) There are several features in the 
data, most notably the interference side-lobe asymmetry, 
not captured by this idealized model. The asymmetry of 
side lobes can be well understood if, instead of assuming 
square- well confinement, wc consider a smooth potential 
well U\ (x) , giving rise to a WKB wave function form for 
the electronic states in the UW [see Fig. 0(B)]. 

As an example, we model the UW by a symmetric 
potential of the form 

U 1 (x)=E Fl \2x/L\ a . (4) 

The limit a — * oo recovers the case of a square quantum 
well while a finite value of a > 2 defines smooth walls. 
The appropriate choice of a should increase with L. As a 
specific example, we consider the case a = 6. We assume 
that the boundaries are soft enough so that one can use a 
WKB approximation for wave function ip n in Eq. (^) . We 
set n F — 100 for the state at the Fermi energy in the UW, 
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FIG. 3: (A) d\M {+) \ 2 /dV, obtained numerically using the 
WKB wave function ^>ioo(a;) and a = 6 in Eq. (^J). The dotted 
curve shows an approximation [Eq. (||)]. (B) Calculated os- 
cillation pattern G oc Vd\M\ 2 /dV using the numerical result 
from Fig. 0(A). vf 1 = vf 2 — v f an< l A^f = 4n/L. 



which corresponds to electron density of w 100 /im _1 0] 
in the most occupied ID channel, with L = 2 fim. The 
theory again gives two separate contributions |M^| 2 to 
| Af | 2 , which are functions of the variables n^, as before. 
In Fig. [j](A), we plot the results of a numerical calcula- 
tion of d\M^\ 2 /dV as a function of «(+). Unlike the 
square-well case, the smooth boundaries lead to a very 
asymmetric oscillation pattern for positive and negative 
re W . Furthermore, the amplitude of the oscillations in 
Fig. p|(A) drops much slower than in the case of infinitely 
steep walls. Results of a numerical calculation of G (V, B) 
for the soft-confinement model, obtained from Eq. (^), 
are displayed in Fig. p|(B). They are very similar, at least 
qualitatively, to Fig. g(B). 

The numerical results can be understood analytically 
as follows. To find M^ + \ we substitute e lk2X for tpk 2 
in Eq. (^), and replace ij) n by the right- moving WKB 
component k{x)~ 1 ^ 2 e lkF ^ x e~ lslyX \ where k(x) = k Fl [l — 
Ui^/E^] 1 / 2 and s(x) = f* dx'[k Fl - k{x')]. Using 
k2 — Q = k Fl — re , we see that if < < k Fl , there 
are two points, x + > and x" < 0, where k(x) — ki — Q 
and the integrand in Eq. (|3j) has a stationary phase. 
We evaluate the integral near these two points using the 
method of steepest descents, and add the results to ob- 
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tain an estimate of M^ + \ In the case of a symmetric 
potential, U\{x) = U\(—x), one finds 



U{(x+) 



■ cos 



(5) 



In Fig. |(A), we plot d\M^\ 2 /dV, obtained using 
Eq. (|5|) as a dotted curve: it is in an excellent agree- 
ment with the full numerical calculation (solid line) for 
large positive k^ + \ The stationary phase approxima- 
tion (SPA) is bad for small values of where the 
conditions for its validity are violated. For the poten- 
tial of Eq. @, 2x+/L « pK^+^/fcp] 1 / for rW < fc P , 
and the phase [n( + ' > x + — s(x + )] in Eq. (||) is equal 
to K^x + a/{a + 1). For negative k^ + \ the equation 
k(x) — hi — Q does not possess a real-valued solution, 
but after deforming the contour of integration [Eq. (||)] 
into the complex plane, we obtain a complex-valued so- 
lution 2x/L « (2|K(+)|/fc F ) 1/a e~ i7r / Q . The argument of 
the cosine in Eq. (^) then has an imaginary part, which 
causes the magnitude of |M^ + ^| 2 to rapidly fall by a fac- 
tor of w e~ 27r /" over each period of oscillation. The 
approximation [Eq. (jjj)] can be used to predict both the 
period of the conductance oscillations and their ampli- 
tude. The period of the oscillations in the case of smooth 



confinement is given by An 
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2tt/(x 4 



"), rather 



than A«W = 2vr/L for the square-well potential. In 
particular, this means that at low k( +%> the period can 
be significantly larger than the value 2ir/L, which it ap- 
proaches for large k^ +s> . The prefactor 1/U[(x + ) in |M| 2 
falls off less rapidly with in the case of a smooth 
potential than the 1/[k' + '] 2 dependence for the square 
well. For Eq. (@), we find 1/U{(x+) cx [kM] 1 /"" 1 , which 
goes only as for large a or 1/VkW for a = 2. 

(Because of the WKB approximation, we cannot recover 
the square-well result by setting a — > oo in our expres- 
sions.) A larger period for small and relatively slow 
fall-off of the amplitude with increasing are both 
qualitatively consistent with the experimental results. 

An analysis using the SPA and WKB approxima- 
tions can also be applied to an interacting electron sys- 
tem in a pair of wires with soft confinement. The 
conductance is determined by Green's function G 12 = 



([*T 2 *,i(o; ) t),*l/i*,'2(!r / ,0)]). The right- and left- 
mover contributions to G\2 can be approximated by the 
respective contributions CJ)? L for an infinite wire, given 
by LL theory M , multiplied by phase factors of the form 



For a pair of coupled, nearly identical 



wires, the leading singularities in C^ L occur at velocities 
v c and v s , corresponding to charge and spin excitations 
that have opposite sign in the two wires, as the symmetric 
modes are not excited in the tunneling process. 

We find that LL theory preserves the key qualitative 
features in Fig. ||(B), such as the asymmetry of the in- 
terference pattern, but brings about additional features 
resulting from the presence of two distinct velocities. In 



particular, the main dispersion peaks split into two lines 
with slopes defined by v c and v s . As a consequence, 
the overlap of the interference side lobes formed along 
the dispersion slopes creates a moire pattern resulting 
in a periodic modulation of the conductance oscillations 
along the voltage axis. The distance between the corre- 
sponding stripes of suppressed conductance running par- 
allel to the field axis is AV^i ow w 2irhv c v s /eL{v c — v s ), 
which is larger than the period of oscillations, AVfast ~ 
4irhv c v s /eL(v c + v s ). Such a beating phenomenon in the 
interference is seen in Fig. |2] and may thus be a direct con- 
sequence of spin-charge separation in one dimension. For 
the 2 /mi wire, from the data in Figs. ||(A),(B), we esti- 
mate AV slow = 1.8 ± 0.4 mV and AV fast = 0.7 ± 0.1 mV, 
and then find v s /v c — 0.67 ± 0.07. This result agrees 
with estimates made previously in a different regime [H. 
A detailed comparison between theory and experiment 
can thus be used to study electron-electron interactions 
as well as the shape of the confining potentials along the 
quantum wires. Another effect predicted by LL theory, 
and present in the data, is a decrease in the value of G 
at low bias, i.e., zero-bias anomaly. In the limit L — > oo, 
the interference side lobes disappear, and what remains 
for the interacting wires are the split dispersion peaks, 
given by the LL theory of Ref. §. A more thorough dis- 
cussion of LL theory for the double-wire system will be 
given elsewhere Q. 
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